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ABSTRACT 


EXPLICIT AND IMPLICIT CALCULATIONS OK TURBULENT 
CAVITY PLOWS WITH AND WITHOUT YAW ANGLE 

Guan-Wei Yen 
Old Dominion University 
Advisor: Dr. Oktay Baysat 

Computations have been perlormed to simulate turbulent supersonic flows past 
three-dimensional deep cavities with and without yaw. Simulation of these self-sustained 
oscillatory flows have been generated through time accurate solutions of the Reynolds 
averaged complete Navier-Stokes equations using two different schemes: (I) MacCormack, 
finite-difference, (2) implicit, upwind, finite-volume schemes. The second scheme, which 
is approximately 30% faster, is found to produce better time accurate results. The 
Reynolds stresses have been modeled, using the Baldwin-Lomax algebraic turbulence 
model with certain modifications. The computational results include instantaneous and 
time averaged flow properties everywhere in the computational domain. Time series 
analyses have been performed for the instantaneous pressure values on the cavity floor. 
The time averaged computational results show good agreement with the experimental 
data along the cavity floor and walls. When the yaw angle is nonzero, there is no 
longer a single length scale (length-to-depth ratio) for the flow, as is the case for zero 
yaw angle flow. The dominant directions and inclinations of the vortices are 
dramatically different for this nonsymmetric flow. The vortex shedding from the cavity 
Into the mainstream flow is captured computationally. This phenomenon, which is due 
to the oscillation of the shear layer, is confirmed by the solutions of both schemes. 
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Chapter I 
BACKGROUND 


I I Objectives 

It i» often Inevitable to hove a cavity on an aerodynamic configuration, ouch as. 
a weapons bay. a landing gear, or similar recessed areas. Carrying weapons internally 
In high speed fighter-bomhe, aircraft ft ,. PB1II) offer, man, benefits, 
ffrmtur maneuverability, long., combat ran,., reduction of weapon ,„ d 
aerothermal heating prcblems. and reduction of radar delation slgna-j. Hove. 

• -It* Ho. a. high ep«d may be a source * llulah ,„ J, 

'‘Vi' 1 m *> “u>« large fiuctuatlone of pressure and volocity. which 

generate strong acoustic ..... This may dam... th. aircraft structure, the 
devices, and affect the avionics on board 

atrtiar reports In th. refs 2) have provided doscriptlon, .f 

city flow.. Many Investigation., both esperlmcnta, ,„d computational, have been 
conducted to study th. fl.wfl.ld , Mld . , wo .„ d ^imonalonal rectangular cavities 
Some examples of .Utmus calculation, of cavity flow, „„„ „ f , j . 5 ^ 

<*lMlc of tola -catch effort arc footed on further understanding c.f the cavity 
"aw phenomanon for deep cavltl... es,^Uu pressure fiue, nation, ^p 

cavities. Furthermore, it is aimed to provide: 

1. valid computational solutions for deep cavity flow, at .up, , „„ lc 
*. qualitative undersundlng of thr..d,m.„.l<m.l lhtou( , h .„ h .^ 

computational graphics, 
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3. wild computational aolution. for flow past . cavity at 

«. comparison of th. solution, fo, cavity n.». from two different num.ri.ol 

schemes (explicit, finite-difference, MacCormuck scheme and implicit, finite- 
volume, upwind scheme), 

6. a study of three-dimensional unsteady flow separation, 

a comparison of time averaged end instantaneous properties of the flow 
solution obtained computational iy with the experimental date. 

A literature sure., on cavity flows pi,.,, „« tw(> The 

Phys'mi of cavity flows ami the,, classification .r, Chapter 2, The „„r„,„, 

equations and the corresponding boundary conditions are detailed In Chapter 3 , The 
solution algorithms and the grid generations are discussed In Chapter 4. The results 
•cussed in Chapter 5, and conclusions and recommendations are given in Chapter 


1.2 Experimental Literature on Cavity Plows 
loshko (21 studied the tlme-averaged effect, of flow, over various cavities |.w 
numbers. Pressure clients various point, on th. cavit, wall. ^ fw 
-ere mealed and friction mmfflcionm wore calculated. Per . deep cavity, . ringl. 
vortex ... men to exist inside. and smaller s~.nd.ry .mrthm. were aim, obmrved a, 
the comers of the cavity. For . slmllm, cavity, ,h. shear lay., w~ ~a„ be 
Inched to the floor of the cavity enclosing vortices on elth„ ride. Roshko conclude 
that the Increased drag due to ,h. eavl„ stemmsd th . , U gn.„on pressure M 

Mm downstream cavity wall as the shear layer impinged on It. ffar.mcheti |„ 
fhrmml on extensl« stud, .f the sound radiated by flow, . v . r mcUngular cavities 
Using a fl., pi.,. »„„ . cvlty of variable l. nph mounrnd In a blow dm™ Ut„„.| 
The Mach number was varied from 0.25 to 1.5. Schllcran ami lnmrf.rom.ler Uchnlqu., 
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wer « used to visualise the acoustic fluctuations and to get quantitative measures of the 
decibel level and directionality of the radiated sound. The boundary layer was changed 
between laminar and turbulent regimes using the flat plate angle of attack and a trip 
wire, and its characteristics were measured using hot wire anemometry. He observed 
that below a minimum cavity length, the shear layer jumped across the cavity without 
impinging on the downstream wail or generating acoustic oscillations. Above this 
length, the wavelength of the oscillations was proportional to the cavity length. When 
the boundary layer was turbulent, several harmonics were observed. However, the peak 
amplitudes observed were considerably less than in the case of the laminar boundary 
layer. 

MacDearmon [61 has presented the results of systematic variations of depth, span, 
and upstream and downstream Up radii on the flow characteristics in a rectangular 
cavity on a flat plate at a Mach numb* J 3.55. Plumbee et al. (7) studied the 
acoustic response of large cavities in flows of Mach number 0.2 to 5.0. Both discrete 
tones and broadband noise generated by the cavity flow were observed. Rapid 
fluctuations occurred in the separated boundary layer, and depending on flow conditions, 
both expansion waves and shocks were observed at the separation point near the cavity. 
The static pressure in the cavity increased with increasing cavity depth. The deep 
cavities would resonate primarily in the depth modes, while longer cavities would show 
lengthwise modes and random buffeting. The most significant drawback of his theory 
was that shear layer turbulence was considered to be the driving mechanism for 
periodic oscillations, contrary to the observations of Karamcheti 1 1 1 and others that the 
oscillations were much stronger when the shear layer was laminar as opposed to 
turbulent. 

Bast [8] observed that deep cavity resonance appeared to result from simultaneous 
doubly -tuned amplification of shear layer unsteadiness by both t'ie shear layer edge tone 
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•nd the cavity enclosure acting ac an acoustic resonator. Three-dimensional effects due 
to variations in the transverse dimension of the cavity did not appear to alter the 
unsteady effects. Rossiter |9| performed wind tunnel tests on subsonic and transonic 
flows ( 0.4 < M < 1.2 ) over shallow cavities. Cavities of L/D >10 were found to 
generate smooth broadband pressure fluctuation spectra, as opposed to cavities at L/D 
of 1 and 2, which generated dominant periodic fluctuations. Higher harmonics of the 
discrete tones were attributed to the distorted wave formes of the oscillations. High 
speed shadowgraph motion pictures showed that the shear layer rolled up into discrete 
vortices, shed periodically from cavity leading edge. 

Spec [101 observed periodic inflow and outflow of air close to the trailing edge, 
and an accompanying lateral displacement of the shear layer. Continuous production 
of vorticity occurred at the leading edge, as opposed to discrete vortex shedding. The 
"captive vortex” or recirculation rone - in the cavity was seen to grow and shrink 
periodically. Root-mean-square pressure fluctuations reached 40% of dynamic pressure, 
which was too large to be explained by the linear acoustic theory. At low Mach 
numbers, the acoustic source exhibited monopole behavior. The radiated sound waves 
were sinusoidal initially, later becoming sawtoothed as they propagated into the far 
field. 

McGregor and White [11] measured the drag of rectangular cavities in supersonic 
and subsonic flow (0.3 < M < 3.0), they also provided schlieren pictures. The pres- 
sure wave generated by the impact of the vortex at the trailing edge went upstream, 
and another vortex was shed. The steady-state drag of the cavity was attributed to 
the impact pressure at the downstream wall. When resonance occurred, a 250% 
increase i drag was observed, which was attributed to the large deflection of the 
shear layer and resultant loss of momentum. Heller, Holmes, and Covert (12) 
conducted wind tunnel tests on cavities with L/D ratios of 4 to 7 over a Mach number 
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range from 0.8 to 3, obtaining spectra of pressure fluctuations inside the cavity and 
under the approaching boundary layer. They also observed stronger spectral peaks 
when the boundary layer was laminar. Pressure fluctuations were seen to be highest 
near the trailing edge, falling offT inversely with distance towards the leading edge. 

Stallings and Wilcox |13| have measured flow past various cavities to obtain 
cavity pressure distributions for a wide range of supersonic Mach number. The test 
Mach number were varied from 15 to 1.86 for cavity depths and widths from 0.6 to 
2.5 inches and cavity lengths from 0.5 to 1.2 inches. These pressure distributions 

together with Schlieren photographs were used to define the critical values of cavity 

length to depth ratio (IVDJ^ that separate deep cavity flows from shallow cavity flows. 
It was further observed that a large improvement in the correlation of measured cavity 
centerline pressure distribution for cavities of various depths were obtained when both 
the cavity width to depth (W/D) ratio and length to depth (L/D) ratio were held 
constant rather than L/D alone. Decreasing cavity width resulted in a reduction in 
(L/D) <r . Three-dimensional effects in the form of large lateral pressure gradients 
occurred on the rear face of the cavities that had closed cavity flow fields. 

Heller and Bliss |141_ conducted wind tunnel tests on cavities of LVD ratios from 

2.3 to 5.5 and over a Mach number range of 0.8 to 2.0. Detailed information on the 

normalized levels of the first three resonant modes in the cavity for a range of cavity 
L/D and free-stream Mach numbers were obtained. It was shown that the cavity 
temperature was close to the stagnation value with recovery factors between 0.8 and 
0.95. Shaw [15] performed wind tunnel tests on a cavity model with variable length 
to depth ratio. The test Mach numbers were 1.6 to 2.86 and Reynolds numbers from 
1.0 to 4.0 million/foot. The model was tested at two angles of attack, two yaw angles 
and two cavity widths. Acoustic data were obtained for almost all combinations of the 
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test parameters. Sound pressure levels as high as 165dB were measured. Reynolds 
number was shown to strongly affect excitation of specific resonant modes. 

1.3 Literature on Prediction Methods for Cavity Flows 

Weiss and Florsheim 116) modeled a low Reynolds . number flow In a deep cavity, 
neglecting convention of vorticity, and obtained results showing steady, double-celled 
recirculating flow inside the cavity. Pan and Acrivos 1 17) used a relaxation technique 
to obtain creeping-flow solutions for rectangular cavities of L/D ratios 0.2 to 4. 
According to their results, the steady flow in a cavity at high Reyn<>)ds number 
consisted of a single core of uniform vorticity, with viscous effects confined to their 
shear layers near the boundaries. For infinitely deep cavities, the viscous and inertial 
forces would be comparable at all Reynolds numbers. 

Mehta and Lavan 118) calculated the flow in a two-dimensional channel with a 
rectangular cavity in the lower wall and the upper wall moving at uniform velocity. 
They solved the Navier-Stokes equations for a laminar incompressible flow. O'Brien (19J 
studied closed streamlines associated with the channel flow over a cavity. The viscous 
Stokes flow in a rectangular cavity with parallel shear flow was calculated by a direct 
finite difference technique, Nallaswamy and Krishnaprasad (20) studied steady cavity 
flows at high Reynolds numbers ( 0 < Re < 50,000). Three fully viscous eddies were 
found inside the cavity, and obtained velocity, temperature, vorticity, and heat flux 
profiles inside the cavity. 

Bilanin and Covert [21] assumed that shear layer instability as well as interaction 
between the shear layer and the cavity trailing edgo were required to sustain discrete 
frequency oscillations. They modeled the shear layer as a vortex sheet excited at the 
leading edge by a periodic pressure pulse, and the pressure fluctuation at the trailing 
edge as an acoustic monopole. They obtained quantitatively correct acoustic mode 
shapes and possible excitation frequencies for shallow cavities for 0.8 < M < 3.0. 



Smith and Shaw (22) developed empirical prediction methode for modal frequencies, 
modal amplitudes, broadband amplitudes, L/D effects, Mach number effects, and 
longitudinal distribution of fluctuating pressures. 

Block (23) studied the noise response of cavities of -varying dimensions at subsonic 
speeds. She included the effect of the L/D ratio in Rossiter's model for oscillation 
frequencies, and the maximum amplitude. Iler formula for the maximum-amplitude 
Mach number was of the form: 

M»(1/K)(L/D)/ {4n)l+A (IVD)b-[(L/D) + 0.514)} (1.1) 

where K is the real part of the wave number of the disturbance traveling down- 
stream, n is the mode number, and A and B are empirical constants. Interactions 
betrv^im the depth and length modes were found. 

Hardin and Mason (24) developed a potential flow model of two-dimensional cavity 
flow in which the shear layer was represented by discrete rectilinear vortices in order 
to predict and explain the broadband noise generation phenomenon. The spectra and 
directivity of the quadrupole noise source determined by their theory were found to 
compare well with observed results for real aircraft. Borland (25) solved the two- 
dimensional Euler equations for the time-dependent inviscid compressible flow over a 
cavity. He modeled the shear layer oscillation and used a piston at the rear bulkhead 
to simulate the mass addition and removal at the trailing edge of the cavity. 

Hankey and Shang (26) analyzed pressure oscillations in an open cavity using the 
unsteady Navier-Stokes equations. With supersonic flow outside the cavity, the Mach 
number inside the cavity was found to be 0.5. The sound inside the cavity propagated 
upstream at about half the freestream velocity. The shear layer was soen to be 
unstable for low frequencies such that 

l/( 4;r) > fK/U (1.2) 

where K Is the boundary layer thickness, and U the freestream velocity. Short cavities 
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who* lengths were lees then 2,tK would not resonate. Above M*28. no Rayleigh 
Instability was found for the shear layer. Peak amplification occurred at half the cut- 
off frequency, so that higher harmonics could exist. The maximum intensity of pressure 
oscillations occurred at about Msl. 

Brandeis 127) studied* the effects of altering the length and aspect ratio of 
rectangular cavities on the development of the shear layer. He used an Interactive 
method which adapted the compressible boundary-layer model for the flow within the 
cavity. His results showed the location of the stagnation points to be sensitive 
Primarily to span variation. Shaw et a! (28) modified Rossiter's formula to improve 
correlation with measured data for cylindrical and rectangular cavities over the range 
0.4 < M < 1.2. Pathasaraty and Cho (29) developed empirical design equations for 

the dimensions, frequencies, and root-mean-square pressure amplitudes of cylindrical 
aerodynamic whistles. 

Baysal and Stallings (30) performed calculations for two dimensional cavities over 
an L/D range of 6,32, and 16 at a Mach number of 1.6, and compared the 
computational and experimental pressure data. An upwind, finite-volume scheme was 
used to solve the complete Navier-Stokes equations. Rizseta (31) presented a numerical 
solution for the unsteady flow over a three-dimensional cavity at a frees tream Mach 
number of 1.6 and Reynolds number of 1.5 million. The self sustained oscillatory 
motion within the cavity was generated numerically by integration of the time-dependent 
compressible three-dimensional Navier-Stokes equations. Comparisons with ««p»rim tn tal 

data were made in terms of the mean static pressure and overall acoustic sound 
pressure levels within the cavity. 

B«)«l. Srlnivasan, .nd Stelll. 1 * (321 ptrfonwd cokul.tlon. thr^lwMloM, 
*"* and shallow cavities of l/D* «. sod 16 for supersonic do,. The, assd tho 
McCormack schomo to sol., tlmadap.nd.nt complete Navior-Stokos 


Comparison* with experimental data were made in terms of SP!U (sound pressure level) 
and time-averaged pressure coefficients. Computational pressure fluctuations wore 
transferred from time domain to frequency domain using the Fourier transformation. 

Suhs and Jordan 1 33 1 used the Chimera scheme to divide the computational 
domain into two overlapping grid regions, then used an implicit algorithm to solve the 
thin-layer Navier-Stokes equations. Numerical computations were performed at Mach 
numbers of 0.74, 0.95 and 1.6, and L/D ratios of 4.5 and 9.9. They have added 
consideration of acoustic suppression devicetsaw tooth fence) into their code and got the 
computational sound pressure level results. 

Baysal and Srinivaaan [34] presented a computational investigation of subsonic 
and transonic flows past 3-D deep and transitional cavities. Computational simulations 
of these self induced oscillatory flows have been generated through time-accurate 
solutions of the Reynolds averaged full Navier-Stokes equations, using the explicit 
MacCormack scheme. The computational results, which are compared with the 
experimental data, include instantaneous and time averaged flow properties everywhere 
in the computational zone. Time series analyses have been performed for the 

instantaneous pressure values on the cavity floor. 

Om [35] conducted a numerical study of cavity-flow phenomenon on a modified 
Boeing 767. A two-dimensional Navier-Stokes code was used to simulate the flow field. 
The code employed the explicit MacCormack scheme. The investigation was aimed at 
examining the unsteadiness of the shear layer and obtaining details of the flowfleld. 
Cavity flow was simulated for two different cavity sizes as well as for two different 
ramp shapes. The computational results indicated that the shear layer stability depends 
very strongly on the shape of the aft ramp, and the effect of cavity size was negligible 
on the shear layer stability. 
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B»y«al and Srinivaaan (361 performed numerical simulations of superaonic turbulent 
flow* by solving the Reynolds-averaged full Navier-Stokes equations by an implicit finite* 
volume method. Several attamples of two dimensional solutions are given to illustrate 
the k-S turbulence model for well and free turbulent shear flows. 


CHAPTER 2 

PHYSICAL DESCRIP" ON 


2.1 Type* of Cavity Plow* 

A city In . flat plate change, the flow field and cue a algnlflcnt Increac 
in the fluid drag. Th. flaw .tractate within a rectangular city I. cmplea. crtkal. 
three-dlmencianal and uaateady. Acceding la th. rati, af th. city length U> depth. 
Uc ca.lt, cn be deified a. dcp. tr.n.itlen.1 and ahallaw 1131 a. Mlceud In 
Plg.2.1. A dcp city la an. fa. which L® I. lac than 10, and crrttlc with L® 
greeter 13 are died .hallow. Wlten L® I. hatwen 10 and 13. It I. cllad a tranal- 

tional cavity. 

par a .hallow cltyf L® > 13 I. the flow U likely «« rctuch on th. bottom 
a the city el the. fane two rclrculattag region. Far a dcp cltylL® < 10). 
a .Ingle vertex U oberved. uaually cmawhat dawn.tre.rn af th. middle af th. city. 
Small ccndary vortic. retattng cunt., t. the main varte. have ben oberved car 
Urn front comer*. In vary dcp critic, a ccnd tmrtc U llkaly ta farm below the 
„r.t ac. forming a double-cllad .tructur. 137). The remnant mad. fmquenclc or dcp 

cavitie* have been observed to change considerably with 1VD. 

The general flow structure of a deep cavity is shown in Pig. 2.2. After the initial 
n pd- af the flow Into the city, a .her layer b formed bat.cn the hlgh-.ped 
external flaw and th. .lower Intarnal flaw. Th. water table rimulaUtm .tudy af th. 
dcp cltle. by Halle and BUc 1141 .haw* that uctcdy matton af the .he layer 
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lead* to mass being pumped In and out at tha cavity trailing edge Thla affect 
generates forward travelling waves In the cavity that reflect from the front bulkhead 
and become rearward traveling waves The shear layer is pulled continuously into and 
then pushed out of the cavity due to the pressure oscillations inside the cavity. When 
the cavity pressure is lower than the freestream pressure, the shear layer is deflected 
downwards pumping mass and momentum into the cavity. This ingested mass Is slowed 
down by various dissipative processes within the cavity, therefore increasing the cavity 
pressure above the freestream pressure. The shear layer is then deflected out of the 
cavity by the excess pressure and mass is pumped out of the cavity with low mo- 
mentum. This oscillatory process extracts additional freestrean. momentum during the 
cycle. A series of sketches which shows the stages of such a cycle are shown in 
Fig.2.3. The choice of a starting point for the cycle is arbitrary, and it is necessary 
to review the entire process to understand completely the conditions at the beginning. 

2.2 Flow in a Deep Cavity 

Cavities with L/D < 10 exhibit open cavity flow structure. In the external flow, 
the wave patterns at the leading and trailing edges of the cavity depend on shear layer 
positions. A supersonic flow past a cavity is shown in Fig. 2.4. High pressure ahead 
of the rear face within the cavity, venting into the low pressure region downstream of 
the front face, cause the shear layer to flow over the cavity. The pressure coefficients 
over the cavity floor are slightly positive and relatively constant with the exception of 
a small adverse pressure gradient ahead of the .wr fr ee, that is associated with the 
shear layer impinging on the outer edge of the rear face. 

In Fig. 2.4(a), a separated shear layer is shown to approach and flow over the 
trailing edge. Thereafter, the shear layer deflects down, thus exposing the trailing edge 
to the free stream, which causes a shock front to occur at the trailing edge. 
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Simultaneously, the downward motion of the sheer layar cautaa a pressure wave to 
appear In the corner region of the trailing-edge bulkhead and the floor. This preaeure 
wavo moves toward the leading edge. The pressure wave also causes the shear layer 
above to bend outward Into the free stream. This outward band causes a shock front 
to propagate, along with the pressure wave, in the upstream direction. 

When the shear layer is above the trailing edge (Pig. 2.4b), there seems to be 
no significant external wave system m a downstream wave is arriving. U must be 
recalled that the shear layer disturbance associated with this wave moves at a subsonic 
speed relative to the external flow. 

At the leading edge (Fig. 2.4c and d), as an upstream wave approaches the leading 
edge, the shear layer is bent downward and there is an expansion wave at the edge 
as well as the upstream traveling compression wave. After the upstream wave is 
reflected to become a downstream traveling wave, the external portion trails away 
since the downstream wave is subsonic relative to the external flow. At this time, 
the shear layer is deflected upward and a leading-edge oblique shock occurs. This 
sequence of occurrences can be compared with observed radiation pattern to arrive at 
Pig. 2.5 for supersonic flow. 

2.3 . ? low in a Shallow Cavity 

The basic flow field in the shallow cavity can be divided into two regions. The 
first half of the flow behaves like the flow over a rearward facing step and the latter 
half is similar to that of forward facing step flow, as illustrated in Pig 2.1a. 

In an example of "Breakaway Separation" (wherein separation occurs on convex 
corners even through there is a favorable pressure gradient), the upstream boundary 
layer is unable negotiate the sharp corner and leaves the wall at the corner, splitting 
into two regions: (1) a slow recirculating separated region near the base, (2) a free 
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shear laysr that ia awapt downstream and reattaches to the cavity floor. At the flow 
meats the 90 degree turn at the front face, it undergoes e rapid expansion. The 
pressure within the separated and recirculation region fails much below the free stream 
pressure. This effect and the pressure of tho side wall enhance the cross flow which 
in turn affects the flow separation characteristics in this region. The separated shear 
layer curves downward and Impinges on the cavity floor. The boundary layer grows 
from the reattachment point up to a certain extent when the adverse pressure gradient 
from the rear face causes separation close to rear face. 

2.4 Unsteadiness and Acoustics of Cavity Plows 
There are two different type of pressure fluctuations that can occur within cavities, 
a random pressure fluctuation and a strong periodic pressure fluctuation. Karamchetl 
(1) found that the periodic pressure fluctuations are accompanied by strong acoustic 
radiation from the cavity. The frequency of these pressure fluctuations was found to 
Increase with airspeed and decrease as the cavity length was increased. In very deep 
cavities (L/D = 1 and 2), the pressure fluctuations are mainly periodic but as the 
cavity depth is decreased, the fluctuations become random in character. For the shallow 
cavity (L/D > 10), the spectrum is smooth and covers a broad band of frequencies 
showing that the pressure fluctuations are random in character. As the depth of the 
cavity is increased, peaks occur in the spectra indicating that periodic pressure 
fluctuations are superimposed upon the random levels. The random component is seen 
in cavities where L/D >4 and the periodic component dominate in cavities where L/D 
< 4 [63). The random component is most intense near the rear wall of the cavity, but 
for very shallow cavities a local region of intense pressure fluctuations occurs where the 
flow reattaches to the floor of the cavity. 
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Bilanin and Covert |21| Analyzed a model of the excitation mechanism suggested 
by Koialter |0|. They treated the pleasure fluctuation at the trailing edge as an 
acoustic monopole As indicated In Chapter I, the mass addition and removal process 
at the cavity trailing edge Is caused by unsteady motion of the shear layer. This 
process produces a piston-like effect at the rear bulkhead (Fig. 2 6), which sets up the 
internal wave structure that forces the shear layer. 

In flow visualization photographs, the shear layer has been observed to roll up 

into vortices that travel rearward and Impinge on the cavity trailing edge, that would 

cause the cavity oscillation mechanism. In supersonic flow experiments, discrete vortices 

are not usually seen, but an amplifying sinusoidal motion of the shear layer traveling 

toward the trailing edge is often evident. Perturbations to the shear layer indicate a 

varying vorticity distribution along its length. The apparent traveling wave motion 

along the shear layer is a result of the dominant amplification of the downstream 

traveling wave in the cavity. Heller and Bliss (14) proposed that the shear layer 

which is subjected to forcing by the cavity internal wave structure would roll up into 

downstream traveling vortices. Thus, the appearance of discrete vortices on the cavity 
♦ 

shear layer is completely cpnsistent with the oscillation mechanism. They conclude that 
"vortex shedding" is a manifestation tf the oscillation process, but it is not essential 
to the underlying mechanism. Fig.2.Y shows the feedback loop responsible for cavity 
oscillations. 
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Fig. 2.5 Typical cavity external radiation patterns 





Fig. 2.6 Pseudopiston 
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Chapter 3 

MATHEMATICAL MODEL 


3.1 Governing Equations 

The governing equations are the three-dimensional, time dependent, compressible, 
complete Navier-Stokes equations in terms of Reynolds averaged variables. These equa- 
tions are written in general coordinates and conservative form [38] as 

^U + ^-(F-Fv) + ^(G*Gv) + ^jH-Hv)*0 . t3.1) 


where U is the vector of conserved variables, 

U * j-[p.f*».pv.pw,pE] T 
and F, G, H are the inviscid (lux vectors, 


F = 1/J 


Vi 


Vs 


“pu, 

puU,+3i5 p 


puu 2 +a^ 2 p 


PuU 3 t8i 4 p 

pvU, +dj£ 1 p 

, g* m 

pvu 2 +a^p 

, H® 1/J 

pviij+a^p 

pwU t +d 3 ^p 


pwUj+djSp 


pwUj +dj$\> 

(E+p)Uj 


(E+p)U 2 

I 

(E +p)U 3 


(3.2) 


and F y , G v , H v are the viscous flux vectors, 


Ov* 1/J 


“o 


“o 




tu 

Fv * 1/J 


&& tkj 



3k^( utki +vt k2 +wt k ,-qk | 


9k^( ut k , +vt k2 +wt kJ -q k ) 


(3.3) 
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£ (•rtkJ+VTkj+w^,-^ jj 

The contravariant velocity components ( U,, U jf U, ) are deflned by 

»d,5 u+d 2 5 , v +a,5 , W( 

u 2 «d J ^ u +d 2 5 J v+ a 3 5 J w> 

Uj-aiS’u+a^’v+d^V 

The stress and heat flu* terms used In the equations above are given by 

* % 

I 3 /*« » 




* 

where k, n, and m are dummy variables. 

** 4, = V4.. 4, = 4 %. 4, - 4rt4z, { . 

Th. WUI mta , E, and the .nr*,. f<j|k>w , 

E » e + l/2( u* + v* + w* ) 

a * CT 
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The perfect gat taw 


P * (y 1) e = p RT (3 8t 

and Sutherland's laminar viscosity law have been used to complete the system of 
equations. For turbulent flow, the total viscosity p, is defined as the sum of molecular 
viscosity (p,) and the turbulent eddy viscosity (p t >. An algebraic turbulent model is used 
to calculate the eddy viscosity. 


3.2 Turbulence Model 

The effect of turbulence is accounted for through the concepts of an eddy viscosity, 
p t , ant eddy conductivity. In the momentum equations, the molecular viscosity, p ( , is 
replaced by an effective viscosity, p ( : 

p # = p + p t = p (l + p/p) (3.9) 

The study of the turbulent flow and adequate modeling of Reynolds stresses (here 

it will be referred to as turbulence modeling) are current subjects of research. They 

have not yet been fully understood, and up to now, we still cannot find a perfect and 

efficient model. The modeling of turbulence is complicated by the fact that several 
♦ 

length scales exist which control the generation, transport and the dissipation of 
turbulent kinetic energy. A simple empirical model that can estimate compressible 
turbulent boundary layers with separated flow has been selected for this study. This 
bas two-layer algebraic turbulence model was proposed by Cebeci (391. It is based on 
the Boussinesq approximation of modeling the general Reynolds stresses by an eddy 
viscosity, similar to molecular viscosity. This model was later modified by Baldwin and 
Lomax |40|. 

The mixing length model employed in the Baldwin-tamax model divides the shear 
layer into an inner and outer region, and is patterned after a method developed for 
attached boundary layers by Celx-d |39|. It is a two layer algebraic eddy viscosity 
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model In which It givtn by 

n»«. • 




y < y 




y < y 

(Twmvir 


(3.10) 


In the inner region, the eddy viscosity is given by the Prandtl-Van Driest formulation, 

P, * (kyD)*|o)| o n) 

where 

D » 1 - exp ( -y*/A* ) 




In the outer region, Baldwin-Lomax model uses the Clauser formulation for the 
outer region, 


P = KC F F (y) 

* mi, KUEB' 7 ' 


(3.12) 


and 


(3.13) 


(3.14) 


P.„ t - min I C -y ^U.JF__ | 

The quantities and are determined from the function: 

F(y) = y|W||l * exp(-y i/A')|. 

In wakes, the exponential term of above equation is set equal to zero. The 
quantity U the maximum value of Piy) that occur, in a profile and y — ia the 
value of y at which It occur.. The Kl.bonolf intermittency correction and U are 

4if 

given by 

P xtss ( y ) = 11 + (3.16) 

U *r = » u *+v* + w* ( u*+vs + wa Ja*^ (3 , 6) 

Th@ Kcond term in is taken to be zero (except in wakes). The values of the 
constants are: 

A. = 26, C_ = 18, C„_ - 0.3, C_ = 0.23, 


P = 0.72, P, = 0.0, C. 
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Tht outer formulation can be used in wakes as well as in attached and separated 
boundary layers. The distribution of vorticity is used to determine length scales so that 
the necessity for finding the outer edge of the boundary layer (or wake) is removed 

The inner and the outer eddy viscosity models are combined to form the turbulent 
eddy viscosity p t in the following manner. First, profiles of p and p # are obtained on 
each coordinate line emanating from the wait. The first point near the wall at which 
P ( exceeds p # is denoted the "cross over point". The turbulent eddy viscosity p t is then 
equal to p t for all points between the wall and the crossover point, and it is equal to 
p # for all points above and including the cross-over point. 

Some modifications have been done to the model for all the points within and 
close to the cavity. These are Degani-Schiff modifications, the multiple-wail modification, 
and the wake modifications. 

3.2.1 Degani-Schiff modifications |41) 

The major difficulty encountered in applying the Baldwin- Lomax turbulence model 
to vortical flows is the evaluation of the length scale, y , and in turn determining 
(*\)«aur boundary- layer profiles in the cross flow separation region. In Fig. 3.1a, a 

general F(y) curve is shown. If there is a strong overlying vortex, then F(y) curve may 
switch to Fig. 3.1b. In addition to a local peak in F(y) in the attached boundary layer 
at y st a, the overlying vortex structure causes a larger peak in F(y) at y = b. The 
choice of the peak at y = b results in a value of F ^ and, in turn, a value of (p^)^ 
which is much too high. Thus, in general, the computed eddy-viscosity coefficient in 
the cross flow separation region behind the primary separation point wiil be too high. 
This will cause the details of the computed flow to be distorted or washed out. 

To eliminate these difficulties we have modified our implementation of the 
turbulence model. At each computational coordinate the code searches outward, sweeping 
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3.2.2 Multiple-wall Modifications |4J 
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3.2.3 Wake model Modifications 

Tho abrupt change in the eddy viscosity coefficient from the boundary layer to 
the wake is approximated using the concept of a relaxation eddy viscosity model to 
represent the different length scales in the problem. 

For the cavity, the eddy viscosity is calculated by 

* P* + <l* u • !»„> II • axp (-x/10 II (3 lg) 

where 

M u : ^e value at the upstream lip, 

M u : Baldwin- Lomax eddy viscosity value, 

(b : instantaneous boundary layer thickness at upstream lip, 
x : streamwise distance from the corner. 

This modification is known as the relaxation model, which has been shown to 
work well for other numerical calculations 142] It accounts for the history effects of 
the fluid. The eddy viscosity (p t ) has been set to zero at all the solid surfaces. 

3.3 Initial and Boundary Conditions 

In computational fluid dynamics, the initial conditions usually correspond to a real 
situation for a transient problem. In practice, initial conditions are obtained from 
experiments, empirical relations, approximate theories, or previous computational results. 
An improper initial guess may result in solution failure. An important requirement for 
the initial condition is that they should be physically as close as possible to the actual 
nature of the flow field. 

A reasonable approach has been to initialize the entire flow field above the cavity 
with the inflow conditions. Within the cavity region, depending on the type of cavity, 
the initial conditions are different. In deep cavities, the flow within the cavity is 
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subsonic which was known from experimental observations- Therefore the velocity inside 
the cavity has been arbitrarily specified as 10% of its free stream value In shallow 
cavities, the shear layer impinges on the cavity floor. As a result, an approximate 
velocity profile close to the inflow velocity profile has been specified within the cavity. 
The pressure and temperature within the cavity has been set to the free stream value. 

Correct boundary -conditions are essential to the success of numerical calculations. 
Six faces require attention in the specification of boundary conditions. On solid surfaces, 
a no-slip boundary condition has been used so that all the velocity vectors vanish In 
the test cases, the wall has been considered to be adiabatic. The pressure at the solid 
surfaces has been obtained by a zero-order extrapolation from the interior point value 
of pressure in a direction normal to the wall. The density is obtained from the state 
equation. 

u = v = w = 0, 6p/6n = 0, 6T/6n = 0. (3.19) 

In the case of supersonic in flow (except in the subsonic portion of the boundary 

♦ 

layer close to the wall), all flow characteristics point from the outside towards the 
inside of the computational domain. Therefore, all elements of the primitive variables 
(vector U in Eq. 3.1) have been specified by a profile generated using the two- 
dimensional compressible boundary-layer equations for perfect-gas flows (the governing 
equations are solved by an iterative three-point implicit finite-difference procedure) (43). 
In the subsonic and transonic inflow, all flow characteristics, except one, point from 
outside towards the inside of the computational domain. Therefore, only the variables 
u, v, w, and T at the inflow are specified. The pressure has been extrapolated from 
the computational domain, to allow for the information to be propagated upstream. 
The pressure within the boundary layer is maintained at the value of boundary layer 
edge. 
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supersonic inflow: P = P , 

subsonic inflow: P (j = P fj (for j above the boundary layer), 

= P # (for j within the boundary layer), 

P # : value of pressure at the boundary layer edge. 

When the outflow field is supersonic, the flow characteristics point from the inside 
of the computational domain to outside. Therefore, ail variables at the outflow can be 
determined from the interior flow solution by a zero-order extrapolation. When the 
outflow field is subsonic, since there is only one incoming characteristic, only one analy- 
tic boundary condition is required. Therefore, we only specify static pressure, and other 
variables are obtained by a zero-order extrapolation, 

bu/6s = 0, 6v/6s = 0, 6w/6s = 0, dT/ds = 0, (3,21) 

dP/ds = 0, (supersonic outflow) 

= K.—’ (subsonic outflow), 

where s indicates the streamwise coordinate. 

The far field refers to the flow field at a distance away from the body which is 
greater than the reference length of the cavity. The boundary conditions are specified 
by zero-order extrapolation from inside * the computational domain for outflow, and as 
free stream conditions for inflow. The pressure values are assumed to be free stream 
at this boundary. 

du/dn * 0, dv/dn * 0, dw/dn * 0, (3.22) 

dT/dn = 0, P = P . 

On the symmetry plane, the following boundary conditions are used: 

du/dn = 0, dv/dn *0.. w * -w, (3.23) 

dT/dn = 0, dP/dn = 0. 



Chapter 4 

NUMERICAL PROCEDURE 


A •thematic of . city model j. shown in Pig. 4 1. It he. Uo„ felllM , from 

oxporimontai observations, that cavity flow, am almost .ymmotric about the cavity 

centerplane (•traamwiae direction). In order to reduce the computational time and 

memory a half span cavity ha> been modeled. One of the modela it shown in Pig. 
4.1b. 

The neat step la the grid generation, which depends upon the flow, the geometric 
configuration, computattonal rnourem. th, mgulrnd accunc, of the solution. the phy.lcal 
dimensions, mul the bounder, layer thlckmss. In the prewnt mull, da, thnmdimen. 
•lonal cartesian grids are used because of tbe geometry. Stretched grids are needed in 
all of the three dlnctlons. so that the entire computational donrnln can be covend with 
a reasoimbl. number of grid points. In which experience large gradient. of the 

no. properties, auch as. don to the mild write, mul shea, layer region, th. grid point, 
are clustered. 

The grids an generated algebnlcally. The cavity is dividml into several sow. 
An exponential stntchlng function given below Is used to cluster the grids in each cone 

y -,-[ (, *‘ 8, -%*k)..)]- ,4,) 

Wh * r *' ttle maximum value of y in the computational zona. The value of k la 

chMon proportional to the clustering matted. The minimum spatial step sics I. of the 
order I0*D, and the maximum cell aspect ratio is around 100. 
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4.1 Explicit, MacCormack, Finite-Difference Scheme (EMFD) 

The MacCormack finite-difference method |44| is a variation of the two-step l.ax- 
Wendroff scheme. This scheme is second order accurate in time and space. This 
explicit method has restrictions on the time step of integration as given by Courant- 
Freidrichs-Lewy 138). It has been used to solve Euler and Navier-Stokes equations for 
numerous practical flow situations, including laminar and turbulent boundary layer 
shock interaction. 

The predictor-corrector explicit algorithm is summarized below in general coor- 
dinates. 

Predictor: 

«r»-u;„-4o; u 

Corrector: 

Au^k- -At t VgOS* ♦ V^k] (4 ' 3) 

Ui.j.k-i[u i(J ,k + Ui.j.k + AU t j'.k] 

where A A,,, A^are forward spatial differences, and V$, V n , Vj are backward differences. 

The coordinate transformation is chosen such that A( ■ At| s s l.o. This 
two step process consists of evaluating spatial derivatives by one-sided differences taken 
in opposite directions on alternate steps. Although the one-sided differences are first- 
order accurate, the combined predictor-corrector step gives the second order accuracy. 

Because of the complexity of the compressible Navier-Stokes equations, it is not 
possible to obtain a closed form stability expression for the scheme. The most success- 
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ful attempt to date has been to linearize the set of differential equations and then 
study the amplification of the Fourier components of the solution by the difference 
method applied to the linearized set. 

For stability, the time step is restricted by the CFL (Courant, Friedrichs, and 
I*wy) condition given by 

CFL-t*».^ (44) 

AlS OiCFL2 
1 + ^Re A 

where X is the maximum eigenvalue of the set of Jacobian matrices, o is the safety 
factor (^0.9), and He is the minimum mesh Reynolds number given by: 

Re A »min (Re*, Re^, Re*) (4 *5) 

Re* 

Re*«£bsdM 

ft 

The linear stability analysis imposes a restriction on the time step in each direction 
(£,t|, and £ ). 


At$< min- 

A* 

i 


At,, < min- 



1 

W|«+ vn,+ wn, 1+ a V( tu+i^mu) 

Atr< min- 

AC 


|uC 1 +vC,+ wC.|+.V(|C!+C^d 


The global time step is determined as At,, 
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For practical computations, especially in the region of large gradients, such as 
•hocks, strong vortices, and separation, the solution exhibits oscillauc- These oscil 
lations sometimes will cause the numerical code to "fair. In ordei to suppress the 
oscillations, we add damping terms known as "artificial viscosity” into the differenced 
governing equation. The fourth order damping term devised by MacCormack and 
Baldwin (45|, is given below: 



where 0 < S < 0.5 and 0 < £ < 0.5. 

The features of the numerical algorithm have been embodied in a solver written 
in 1*0 It i’KAN 200 language using 64-bit arithmetic specifically for the CYBER- 205 
computer 1461- The vectorization executes in | planes. Approximately 3.25 million 
words of memory is required for grid points. The mean data processing rate is 5x10 « 

cpus/time step/grid point on the CYBER-205 computer of NASA Langley Research Cen 
ter 


4.2 Implicit Upwind, Finite-Volume Scheme (IUFV) 

The finite volume discretization is another means of developing difference 
approximations to partial differential equations. It is in this form that the equations 
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express retes of change in a finite volume in terms of the resultant fluxes through the 
boundary of the volume. These equations provide average values of the solution over 
the volume. The advantage of the integral form is that it remains valid in the presence 
of discontinuities in the flow, such as a shock or a separation. The integral form 
(stationary control volume) of the governing equations is 


^/// f U.dV + // t F.n<fc«o . 


(4.8) 


where p = l(P*P v ), (G-G ¥ ), (H-H)| and n«n t f+ n„n + «;C ls * un!t normal vector 
pointing outward from the surface S, bounding the volume V. We can use the 
divergence theorem to change the surface integral to a volume integral 

au - (4 - 9) 

//4?i.dV + ///V.FdV«0 . 

o X 


Assuming continuity of the integrand in the above equation, the differential conservation 
law form of the governing equations can be written as 




(4.10) 


For generality, the governing equations are put into nondimensiona) form. We use the 
cavity depth as the reference length L 
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The finite volume method is adopted to handle arbitrary configurations since the 
method can more easily treat complicated grid structure than finite difference formula- 
tions. We can discretize the flow field into a set of ordered hexagonal ceils (volumes) 
and apply governing equations to each volume directly. Thus, a semi-discrete represen- 
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tatlon of gowning equation can be written at 


[£ <u>v l 


(4.12) 


or, in more detail 


(|SU\ (P* F vl4 ! ik*(F*F»)i.Lik iO-G,) ljt t k .(O.G,) tj .J k 
Idl/Uk A$ At) 


(4.13) 


At) 

( H * H„)y, k *I*( H ♦ H* )i,i k .i 


AC 


I « 0 


where, for convenience 

AC»Ci«l-C^«l * (4.14) 

An*na-ni.» =i , 

2 2 

AC»C»4-Ci-* = i . 

2 2 

and by virtue of the integral representation, 

//**•'•<••>«]*. • <4151 

U is regarded as a cell average value rather than a pointwise value at the cell center. 

For solving the Navler-3tokes equations, thr theory of characteristics is crucial 
in determining the directions of the signal propagation, the information gained from 
characteristic theory. Recently, much interest has been generated in the development 
of numerical m et h ods that model the underlying physics, as dictated by cha ra ct e r is t ic 
theory, at each point. Some of the methods include flux-vector splitting, and flux- 
difference splitting. These methods can be classified as upwind methods and have the 
advantage of being naturally dissipative. The advantages of upwind methods over 
central difference techniques are 1: enhance stability properties, and 2: eliminate 
explicitly artificial viscosity or damping terms. 
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Since the pioneering work of Stogor and Warming! 48) in the early 1980* , the 
upwind method has grown rapidly. The achemee ueed In here are based on the method 
developed by Beam and Warming (49,50i. The scheme for solving the compressible 
Navier-Stokes equations belongs to the same class of AOI methods developed by 
McDonald and Briley (64). 

After Bulor implicit time integration and linearization, the governing equations 
can be expressed in difference form as: 



=jd 4 (F-F^G.G v )+a<{H-Hj]“ -R“ 


These equations are very difficult to solve due to the large banded matrix. The left 
hand side requires an inversion of a very large matrix. There are a number of 
approximate factorizations to split the resulting large banded block matrix equation into 
a sequence of easily solvable equations. We use a three factor block tridiagonal scheme 
that is fully vectorizable in which the implicit operator is factored such that eac h 
directional factor contains the Jacobians associated with that direction only. The three 
factor scheme can be written as 



Apply spatially split approximate factorization 



(4.18) 
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The solution only requires block tri-diagonal inversions end is followed by the update 
step U**> = U» + AU. 

The essence of dux difference splitting techniques 151) is the solution of local 
Kiemann problems stemming from the consideration of piecewise uniform states between 
cell interfaces and an initial data line. Roe |52) proposed a method of exploiting the 
fact that the Riemann solution for any set of linear conservation laws can be easily 
computed. Roe's idea was to obtain an exact solution of the Riemann problem for the 
following linear hyperbolic system, rather than using an approximate solution to the 
exact equation. This system, is approximated by 

|^+[A]||U L ,U t )^=0 . (419) 


For flux terms, a monotone upstream centered scheme for conservation laws(MUSCL) 
approach is used. For example, in the £ direction, the spatial derivatives are written 
conservatively as a flux balance across a cell as: 



where the interface flux i> constructed as. 

F 4-jf p ( q ) +p d q )•! * Ku V)]^ . 


(4.20) 


(4.21) 


and q denotes state variables on cell interfaces determinated from upwind biased 
interpolations of the primitive variables. 


Ri*I“qi+| 1 *k)Z. +•( 1+k) A+J j| , 

Qnl*Qi *| l *«)£♦+( 1 +k)2u] | 


(4.22) 


L- ^63 d -]L 
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where 

(/\ ) s & - q , (A ) s q - q 

♦ • ”«i “ * .i " i 

and the elope limiting ia uted to maintain monotonicity 

A.sfiw[o,ai){A<agnA.,flAjBnA«)]«nA, , 
A.eiw[o,nm(jA.sgnA.,|lA,9gnA.)]sgnA. . 
6 J3-k) 

p cmn ♦ 


(4.23) 


where 


central difference 
flux-difference split 


rewrite the Jacobian metric 


A “^ 3 U = T ( a+ + a ~)t “ 1 , 


(4.24) 


where 


a * s a±JaJ 

2 


are diagonal matrices formed from the eigenvalues of A, i.e. 
Aadag(X t . Xj, Xt.X4.X3), 


Xi “Xt’ 


.X..5JS1. 


v,.!?!it!lsJ . 


(4.25) 


(4.26) 


The contravai t velocity norma! to tha cell interface U 


Bi + + . (4 27) 

and the * Indicates a Roe averaged variable! 63 1 

In order to develop the aolutlon using considerably less computer time, each 
direction of the spatial factor ts approximated separately with a diagonal inversion 
separately, such as 

ijr * t ]’ 4U ‘ * T [iir * * A * ♦ *» v ] tW « » 


both side are then multiplied by T>, and hence the I-sweep then becomes: 


- — + 9* a + 3»a 
JAt * S 


(t*W )«.t*‘r. 


<4.2&) 


Due to the repeated eigenvalues (X, = X t - ^ t )* only three scalar tridiagonal LU 
decompositions are required for each line. The tridiagonal matrix can be written as 


•A {m,.i.u 14 )(tW),.,+ [,i T A*(M,a.u 1 )-.' , (M,^u,)J( T l 4u) 1 (4 
+A (MuJ.U w J(|t‘ , AU*) w «-T^R 1 . 


The metric terms M and state variables U are evaluated at cell centers. After the 
initialisation and development of the flow, the diagonal inversions are turned off and 
more accurate block inversions are turned on. 

The Van Leer type flux-vector splitting method (54,55) distinguishes between the 
influence of forward atid backward moving particles in this model, the interaction 
between panicles is done through mixing of a pseudo-particle going in and out of each 
cell according to a given velocity distribution. 
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Since the flux vector* ere homogeneous function* of degree one In U, they ere 
expressed in terms of their .Jacobian matrices Using a similarity transformation allows 
Jacobian matrices to be written as 166, 67, 301 

(4.311 

f-IaIu-Itat'Iu . 

The matrix is a diagonal matrix composed of the eigenvalues of A and is given by 


A 


X, 0 0 0 0 

0 0 0 0 

0 0 X 3 0 0 

0 0 0 X 4 0 

0 0 0 0 


(4.32) 


where 

^i.a.j*U«^u+^v+5«w , (4.33) 

X^U+lnJa , 

Xj=U.|V*|a . 

The eigenvalues can be decomposed into nonnegative and nonpositive components 

MXT+XT . W.34) 

where 

. ..MM ,4 38 > 

2 

♦ ♦ 

Similarly, the eigenvalue matrix can be decomposed into a«a + a , where a is 
made up of nonegative contributions of a and a is constructed of the nonpositive 
contributions of A . This splitting of the eigenvalue matrix, combined with equation 
(4.7) allows the flux vector to be rewritten as: 

F-t(a%a‘)t 4 U-(a 4 +a)u«F%F' . 


(4.36) 




The flux vector hat three distinct eigenvalues and can therefore be rewritten at a turn 
of three tubvectors, at following: 

F f 

The forward and backward flux vector* K* and F* are formed by inaerting A , = A' 
and A ( = A". retpectlvely. 

For the purpose of splitting the flux vector, approximate factorization it used to 
solve the vector flux split equation, which in delta form, it the following: 


l+AtV 4 Al+AtA t Ai-Atdj(JiLj| x I+AtvX+AtA n Bl-Ata*(!jL 
4 1 ♦ A f ? {d ♦ A t A ? cl - A taj (, j| ( AU ) = -A t{ RHS ) 


( 4 . 37 ) 


KHS«A(F +A$F +A,,G ♦A^O +A;H +A{H • „ 

The split flux vector in the above equations is implemented as the flux across a 
cell, corresponding to MUSCL type differencing: 


V<f\.A*F‘ «[F*( u )+F (|U*)]i*i - [F*(u)+Fdu*)]i.L . 


( 4 . 39 ) 


The notation F«(U*) tl/J denotes the forward flux evaluated using the metric terms at 
the cell interface (1+ 1) and the conserved state variables on the upwind side of the 
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interface, obtained by r, fully upwind second-order etate variable interpolation: 

* 15U, • 0.5U ,, (4 40) 

U * 1.6U • 0.6U 

M • d Ml m3 

All the computational cells are advanced at equal time steps (time-accurate calculations 
require global time stepping) corresponding to a speeded CFL number given by: 


CFLs At( ! U !+| V |+| W |+a[| gnd($ |+| grad(q)|+| grad(0|]} 


(4.41) 


After that, the conserved variables are updated 

= U» + Au (4.42) 

The features of the numerical algorithm have been embodied in a solver written 
using FORTRAN 77 with 64-bit arithmetic specifically for the Cray-2 computer. 
Approximately 10 million words of memory are required for 150,000 grid points. The 
mean data processing rate is 6x10* cpnsftime step/grid point on Cray-2 of NASA tangley 
Research Center. 


' 4.3 Computational Flow Visualization 

The visualization of, the solution through graphics is an important area in 
computational fluid dynamics. This is essential since the amount of data generated by 
solving the equations is in the order of millions. Details of flow visualization have 
been addressed in (581. The PLOT3D |591 is a popular graphics application program 
to create and interactively view flow field solutions. It is an application program in 
that it deals specially with computational fluid dynamics grid and data. In the current 
research, this program has been used extensively to represent the large amounts of data 
pictorially. The contours of pressure, density, temperature, and Mach number have been 
plotted. 
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4.4 Time Series Analysis 

In the current study two characteristic times have been defined for convenience. 
The first characteristic time, t j)t is defined as the time required for a fluid particle 
traveling at free stream velocity to traverse the length of the cavity. This parameter 
is useful in determining the computational time required to get past the non- physical 
initialization of the flow in the computational domain. 

t LAJ 

Cl 

But the acoustic waves propagate within the cavity at local sound speed, so we can 
define a second characteristic time as 


t = L/a = Mt 

«* cl 

where a* is the local sound speed. For reasons of economics, large time steps are . 
desirable for the computations. However, if this time step is larger than t^, then we 
cannot capture the pressure fluctuations. Beside the numerical stability, therefore, there 
is another restriction to the numerical time step. We express that using the frequency 


f > f 

CFO pr*uuf* fluctuntioo 


and from the definition f = Iff, 


^ ^CFD ^ (luctuiVtuo 

= > *cro < roln( V > 


A lot of time dependent data (pressure history) is obtained in the time domain. 
A widely used method is to transform this data from the time domain to the frequency 
domain using Fourier transforms. Compared with other methods, the frequency domain 
representation is much easier to understand and contains the information of harmonic 
resonance frequency directly. 

The definition of Fourier transformation |61) is given as following: 

F( w)»/ T f(l)«pl-iwt)dt . (4.43) 
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There are two fluctuation mechanisms that occur in a cavity; a strong periodic 
fluctuation, and a random fluctuation. Since we cannot predict the period of a random 
fluctuation, we rewrite equation (4.43) as 

F(w)«y fQ0exp(-iwt)dt . (4.44) 

The acoustic pressure distributions are represented as overall sound pressure level, 
SPL, in decibels which is defined us 

I -2 \ (4.45) 

SPL ( t ) * 10 log I J . 

where the value of reference pressure, P, is 2x10* Pascals. Then time distributions 
of SPL are transformed into frequency distributions using Fourier transformations for 
cyclic but non-periodic oscillations (601, 

SPL(f)*^ SPL (t)expf-2itfti)dt . 

(4.46) 

where f is in Hertz (Hz). 

Computed sound pressure levels are corrected by subtracting 

(447) 

to account for the disparity between the computational band width increment, and the 
experimental band width increment. 






Chapter 6 

RESULTS AND DISCUSSION 

The computational test cases have been chosen in order to compare the flow 
structures with different flow conditions (with and without yaw angle), and the 
computational results of different schemes. Four cases have been modeled and analyzed. 

Case 1: Mach=2.16, L/D=3 0, yaw=0», EMFD Scheme; 

Case 2: Mach = 2.16, UD=3.0, yaw =45-. EMFD Scheme; 

Case 3: Mach=165, L/D=6.7, yaw=0>, EMFD Scheme; 

Case 4: Mach = 1.65, L/D=6.7, yaw=0- IUFV Scheme. 

One of the reasons for choosing these test cases has been the availability of 
experimental data. Also, it is intended to demonstrate the capability of the numerical 
schemes to model the complicated flow features in different flow regions. The flow 
conditions and cavity specifications chosen for the four test cases are given in Table 5.1. 
The grid sizes used to represent a field size enclosing each cavity are shown in Table 
5.2. The field dimensions and cavity dimensions are also shown in the X, Y, and Z 
directions for the test cases. 

Since the flows being simulated are unsteady, global time stepping has been used 
in order to maintain the time accuracy. In order to stabilize .the initial numerical 
transients, the time step has been gradually increased by increasing the Courant 
number form 0.01 to 0.7. Shown in Table 6.3 are the characteristic time and the 
eo mpu tet innal schemes, the number of characteristic times over which each of the cases 
has been run, and the amount of CPU (Central Processing Unit) hours used for each 
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structures. For clarity, only ovary other velocity point In the etreamwiee direction hae 
been plotted in all the velocity vector plote. 

In order to demonetrate the unsteady and three-dimensional behavior of the flow, 
the velocity vectors within the cavity at different spanwise planes have been plotted in 
Pig.5.3 at two instants of time <t t| * 40.4 and t f) * 55.76). These velocity vectors are 
presented at four streamwise planes, moving from the plane of symmetry towards this 
sidewall, (Z/W — 0.1627, 0.2728, 0.3672, and 0.4415). It can be seen that the flow 
structure changes from one instant to the next. In addition to the variation with time, 
the spanwise variation is evident The vortex structures near the front and rear faces 
appear to be warped as the side wall is approached. 

The instantaneous (t^ — 55.76) streamwise density contours at four planes (Z/W 
= 0.1527, 0.2728, 0.3672, and 0.4415) are displayed in Fig 6.4a. The shear layer is 
evident from these contours. The value of density at the front face region is low and 
progressively increases to a high value towards the rear free region. The reason for 
the increase in density at the rear free region is the interaction of the shear layer 
with the roar free and the compression of the fluid. This compression causes a free 
stream shock wave close to the rear face. Also, from Mach number contours, we can 
see the shear layer deflects out into the rear face region, and causes a shock wave in 
the freestream. The static temperature within the cavity is slightly higher than the 
freestream value, and the density in the cavity changes from a low value at the wall 
and increases progressively towards the cavity opening until it reaches the freestream 
value. The density decreases in magnitude from the plane of symmetry towards the 
sidewall. The streamwise pressure and Mach number contours are shown in Fig. 5.4b 
and 6.4c, respectively. The same trends observed in the density contours are seen in 
the pressure contours. The pressure contours exhibit shock structure weakening in 
magnitude from the plane of symmetry towards the sidewall. 
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The Instantaneous streamline plots generated with "numerical dye injection" at 
different vertical planes, stretching from the floor vertically up to the shear layer region 
are displayed in Fig. 5.5a. The three-dimensional view of the instantaneous streamline 
within the cavity are depicted here. The flow inside the cavity consists of a large 
counterclockwise routing vortex, and a smaller clockwise routing vortex. We can get 
the deUils of the flow structure from Fig, 5.5a and Fig. 5.5b. 

The instantaneous velocity vectors at three cross sectional planes (X/L = 0.0679, 
0.6852, 0.9571) and two insUnU of time (t t) = 40.4, 55.76) are shown in Fig. 5.6. It 
is seen that the flow structure is not only different at various cross sections but also 
varying with time. The cross flow velocity vectors (at t f( * 65.76) at X/L = 0.0679, 
show the flow on the SFP to be towards the cavity- forming one bigger clockwise 
routing vortex, and a smaller counter routing vortex. The air flow is from the side 
flat plate into cavity. The velocity vectors at X/L = 0.6852, and t , = 66.76 show the 
fluid is being pumped out from the cavity towards the side flat plate. 

The instanUneous limiting streamline patterns on the cavity floor are displayed 
in Fig. 5.7a. Since the limiting streamline ploto do not indicato the direction of the 

stresses, two components of the shear stress contours on the cavity floor ( r and r ) 

w st yt 

are plotted (Fig. 5.7b), and shear stress vectors at the same instant of time as Fig. 
6.7c. The main flow interacU at the rear face and reatUches on the floor close to the 
rear face (Fig. 5.7a). The reattached flow moves in the direction of the front face. 
Towards the front face of the cavity, the limiting streamlines converge into a line of 
separation. There are several nodes and saddle points. This is an instantaneous open 
type separation. In addition to the main separation, there is secondary closed 
separation on front face. This separated flow is seen to reatUch on the floor very close 
to the front face. In addition to the streamwise reattochment, a cross flow re- 

attachment of the flow is seen close to the sidewall. Close to the front face the 


»tre»» vectors ere emell. Since the 


velocity components are small, therefore, the shear 
shear layer interacts with the rear face, the magnitude of the velocity components are 
larger and therefore the shear stress vectors are larger. Due to the fluid particles are 
drawn towards the sidewall in the rear region of cavity, the low pressure region toward 
the sidewall. The shear stress contours are shown in Pig. 5.7b. 

The mean pressure coefficient distribution along the cavity centerline and rear 
face, obtained computationally and experimentally (13), are shown in Pigs. 5.8a through 
5.8f. The specific locations are: (a) flat plate ahead of the cavity (PPI), (b) flat plate 
downstream of the cavity (PP2), (c) front face of the cavity (PF), <d) cavity floor (P), 
(e) rear face in vertical plane (RF), and (fl rear face in horizontal direction (PP2). In 
general, good agreement is seen between the computed and experimental data. The 
computed values of C f on the rear face are lower than the experimental data. Also, 
the computed values of C f on the cavity floor are slightly lower. 

Discrepancies in the computational results can be attributed to the following 
reasons: 

(1). turbulence model, (2). coarseness of the grids, (3). numerical truncation 

error, (4). explicit addition of artificial viscosity. 

Shown in Pig. 5.9 are the frequency spectra of the overall sound pressure level 
at two locations along the cavity centerline, where data has been recorded. The pickup 
points are located in the cavity floor X/L = 0.4, and 0.8. Prom experimental 
observations |I5J, it is known that there are two resonant frequencies. A time series 
analysis in the frequency domain has been performed for histories of instantaneous 
pressure. The pressure data have been converted into sound pressure level (SPL) in 
decibels <dB). The frequency domain has been plotted up to 1.2 megaHerta and 
compared with the experimental data. Basically, the computations predict the resonant 
frequencies, but cannot predict the SPL magnitudes accurately. The discrepancies in 
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the sound pressure level may be attributed to the relatively short periods of data 
collection. 

6.2 Case 2 : Mach * 2.16, L/D * 3, yaw * 460, EMFD Scheme 

Because of the yaw angle, this is not a symmetric flow anymore. So we need 
to calculate the full cavity. In order to save the time for developing the flow, we used 
the primitive variables output file of Case 1 as the initial conditions of this case. We 
can envision the cavity flow with yaw to be a mixed type of two different 1,/D ratios; 
L/D=3 in the streamwise direction and W/D=5 in the crossflow direction. 

The instantaneous velocity vectors (tc = 28.9) of longitudinal planes at four 
different spanwise locations, (Z/W = 0.0456, 0.6250, 0.8301, and 0.9783) are shown in 
Fig.5.10. From this figure, we can see a larger vortex and a smaller vortex within the 
cavity, which is a basic deep cavity structure. The instantaneous density contours 
(tc=28.9) at the same four spanwise locations are shown in Fig.5.11. The organized 
behavior of the shear layer is clearly seen as it bridges the cavity opening. The shear 
layer interacts at the rear face and causes a reversed flow within the cavity, resulting 
in a bigger vortex. At this instant, the density contours show that mass is being 
expelled at the cavity rear face. Due to the expulsion of mass, the shear layer is 
deflected up. There is a large region of separation on the rear flat plate (FP2). This 
is caused by the expansion around the sharp corner at the rear face. The variation 
of properties it. the spanwise direction of Case 2 shows similarities to the typical 
structure of a deep cavity at yaw angle equal to zero. 

Presented in Fig.5.12, are the instantaneous (t c| = 28.9) cross flow velocity vector 
plots at four axial locations (X/D= 0.2874, 0.8676, 2.5165, and 2.9391). At the X/D = 
0.2874 location, the flow within the cavity is going in the reverse direction to the 
freestream, and form a smaller counterclockwise vortex close to the front side wall 


65 


(FSW, sec Pig, 4.1«). In the rear .Id. w.ll, the flu w dip. into the cavity, and deflect, 
out in the reor tide of the cavity. The drn.ity contour, of the cro*. fluw are shown 
in Pig. 6.13. The shear layer i. very .table near the front of the c.vity, but there are 
a lot of fluctuations In the rear part c f the cavity. We still can aee the shear layer 
bridge the cavity. Pressure and Mach number contours are shown in Pig. 6.) 4 and 

Pig. 6.16. Both th.> longitudinal and the cross flow direction figures indicate the 
characteristic, of a deep cavity. 

The instantaneous (t a = 28 9) particle trace, limited to four .treamwise planes 
(Z/W = 0.0466, 0.321, 0.86, and 0 3783) ere displayed in Pig. 8.16. The formation of 
a larger vortex and two counter rotating vortices are evident from these figures. The 
particle traces limited to cross flow plane, at five different location. (X/D = 0.2874, 
0.8676, 1.7998, 2.6165, and 2.9391) are shown in Pig 5.17. The particle traces at four 
different horizontal planes within the cavity (Y/D = 0.0276. 0.2323, 0.5, and 0.8936) are 
shown in Pig. 6.18. We can observe that shear layer brings in air mass and 
momentum. The particle traces are shown in Pig. 5.19. The .hear layer first impinges 
on the rear face, then it deflects toward the rear side wall. It bends downward 
towards the cavity floor! A clockwue rotating vortex is formed wUhin the cavity. 
Prom the top view, clockwise vortex is the main vortex direction. That is the most 
distinctive difference from the zero yaw angle cavity flow. In Case 1, the main vortex 
is predominantly in the vertical plane of the cavity. But In Caw 2, the predominant 
plane of the main vortex inclines at a certain degree with the vertical plane. It is not 
only rotating in the streamwise direction but also in the horizontal plane. 

The shear strew vector, and contour, on the cavity floor are shown in Pig. 5.20a 

•nd Pig. 5.20li. A* expected, the main concentration of the .hear stress is at the 
downstream corner of the cavity. 
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6.5 Cases S «nd 4 ; Machs 1.66, IVD=6.7, EMFD and IUFV Scheme* 

In an attempt to compare the solutions of the EMFD and IUFV scheme*, Case 
a and Case 4 are discussed together. The major similar features between these cases 
are as follows: I. the physical domain and the gvld which covers it; 2. the governing 
equations of tlte flow; 3. the turbulence model; 4. initial and boundary conditions; 6. the 
physical elapsed times (up to two significant digits). The major nonsimilar features 
between these cases are the following: 1. the solution and discretisation algorithms; 2 . 
the manner in which the physical domain is mapped on to the computational domain, 
he., the evaluation of the coordinate transformation metrics and Jacobians; 3. the 
numerical dissipation, i.e., the artificial dissipation of EMFD versus the natural 
dissipation of IUFV; 4. although the total elapsed times are similar, the computational 
time steps which add up to the elapsed time are totally different; 6. 32-bit arithmetic 
(half-precision) of the CYBER-205 computer is used for the EMFD scheme, but 64-bit 
arithmetic of the CRAY-2 computer is used for the IUFV scheme. 

A preliminary numerical study was conducted to determine the effects of eddy 
viscosity. First, a case was run without any eddy viscosity, i.e. as a laminar flow 
Then, the eddy viscosity was included in the computations only on the horizontal walls. 
Both of the test cases, predict the shear layer impinging on the cavity which is 
incorrect for open cavity flow. When the eddy viscosity is computed as mo dified as 
described in chapter 3, then the correct predictions are obtained. 

The instantaneous! t t) = 3.253 ) velocity vectors at two spanwise planes 
(Z/W=0.788l, 0.9764) of both cases are shown in Fig.5.21. In both cases, we can see 
a clockwise rotating vortex which dominates the cavity. The shear layer is clearly seen 
to bridge the cavity in both cases. 

The density contours! t <( =3.253 ) of both cases in various spanwise planes are 
displayed In Fig.5.22. The density contours show that mass is being expelled at the 
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rt«r portion of tho cavity. Due to tho expulsion of man, the »ha«r layer la deflected 
up. Wo can mo a shock generated at tho rear face, and there Is separation on the 
roar flat plate (FP2). 

Tho instantaneous! t <( = 3.253 ) Mach number contours of both coms in various 
spanwlM pianos are shown In Fig 6.23. When the shear layer is deflected up in the 
rear portion of the cavity, we can Me the vortices generated in the shear layer. The 
comparison of instantaneous streamlines! ^=3.253 ) in various planes of both coms 
are shown In Fig.6.24. We can Me a big vortex within the cavity, and the shear layer 
Is deflected out in the rear portion of the cavity. The instantaneous limiting 
streamlines! t cl - 3.253 ) on the cavity floor of both cases are shown in Fig.5.25. We 
can see the flow reattachment in the rear portion of the cavity. Then the flow to in 
the upstream direction, and finally, it Mparates on the front portion of the cavity. 
Also, we can see the three-dimensional separation lines in both cases. In Com 3, there 
is a very strong Mparation generated in the middle portion of the cavity. 

From above figures we can understand that both schemes produce approximately 
similar flow structures in the cavity. 

The instantaneous! t rt = 12.875 ) velocity vectors of longitudinal planes at two 
spanwtee locations! Z/W= 0.4857, and 0.7881 ) are displayed in Fig.5.26. Note that 
theM figures are for t„ = 12.876. which is 9.6 characteristic times after the instant of 
Fig.5.21. The general flow structures are similar. A single vortex dominates the 
cavity. The shear layer interacts with the rear face and it deflects up. The density 
contours! t rt « 12.875 ) of the same spanwtee locations as Fig.5.26 are shown in Fig.6.27. 

When the shear layer is deflected up, shock waves are generated in the rear portion 
of the cavity. 

The Instantaneous !t u ■ 12.676) total pressure contours of longitudinal at 

various spanwtee locations (Z/W=0 , 0.3486, 0.6027, 0.7327) af both coms are shown In 
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- Fig. 6.28. The vortex cells ere seen to exist in the shear layer. The discontinuous 

wrte* is distributed from the lip of front faee to the rear face. As, we also can see 
: from the density contours at two different instants, that is t j( =3 253 (Flg.6.22) and 

J. t ,*12.875 (Fig 5.27), the vortex cells change their positions. Somo vortices are above 

^ the ehear layer towards the rear of tho cavity. So, we can conclude that there Is 

- : vorto * shedding phenomena existing in the shear layer In both cases. That is, It dips 

Sf int0 the cavity and brings mass with It. Then It bridges the cavity for a few instants 

--•j ^th no mass entrainment between the internal and external flows. Finally, when the 

'J ,hear , *y« r deflecto up to open the cavity, vortices are shed from the cavity to outside. 

This phenomenon is confirmed by the solutions of both schemes as evidenced with the 
hotel pressure contours and density contours. Hence it is concluded that the vortex 
shedding phenomenon captured here computationally is a physical one. 
i, The averaged pressure coefficients from both solutions are shown in Fig.6.29. 

i 

The plotted Cp distributions are for the cavity floor centerline, the centerline on the 
-«• rear fac ®> span wise distributions on the cavity floor and the rear face, and the 

longitudinal direction on the side wall. Generally tney show good agreement with the 
experimental date. The predictions of the IUFV method are slightly better than that* 
of the GMFD method. 

; The maximum stable time step for the IUFV code (7.48x10-7) is approximately 1.7 

— '■ times that of the EMFD code (4.39x107). This comparison is for the solutions advanced 

to the same characteristic time. However, each time step computation with the IUFV 
= • code significantly more computer time than the GMFD code. Therefore, a frett - cr 

comparison can be made as follows. The GMFD code needs 4 mega words computer 
memory, and 6.68 CPU hours for each characteristic time. The IUFV code needs 
approximately 10 megawords computer memory, and 4.31 CPU hours for »pi>h 
characteristic time. That is, the IUFV scheme is approximately 30% fester but it 
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requires 2.6 times rr.ore computer memory. Note that the maximum time atep (KM) 
la at laaat one order of magnitude imaller than the characteristic timea t #( (l.fixlM 
and (S.lxlCH). Thia la necesaary to capture the phyaical fluctuations with the 
numerical time accurate solution advancing. 
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Table 5.1. Test conditions for die four cases 


Mach (M) 

T in K 

Re/meter 

number 

static 

in million 

2.16 

166.1 

6.56 

2.16 

166.1 

6.56 

1.65 

210.3 

6.56 

1.65 

210.3 

6.56 


11425 


11425 


Depth 
in m 

Width 

inm 

Length 
in m 

.0127 

.0636 

.0381 

.0127 

.0635 

.0381 

.1108 

.0727 

.7458 

.1108 

.0727 

.7458 


Table 5.2. Cavity dimensions and computational grid 


Case 
No. X 


Grid Size 


Held Dimensions 


Cavity Dimensions 


Block 


Z X 


1 97 60 36 .1076 0587 .0512 0127 .03175 .( 


60 71 .10761 0587 .1024 .0127 .0635 0381 


5 1 120 1 63 28 1 -1076 I .3753 .1163 .7458 .1108 .0727 


120 27 


20 9 


1.189 .3753 
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Fig. 5.24 Instantaneous streamlines at 
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Fig.5.29. Mean pressure coefficient distribution on the cavity along the (a) centerpline of cavity 
, (b) centerline of rear face, (c) crossflow line at Y/D ■ 024 on the rear face, (d) crossflow 
at Y/D ■ 0.68 on the rear face, (e)crossflow line at Y/D * 0.94 on the rear face, (f) crossflow 
at X/L » 0.204 on ih?. floor, (g) crossflow line at X/L * 0.477 on the floor, (h) crossflow 
at X/L » 0.75 on the floor, (i) crossflow tine at X/L ■ 0.885 on the floor, 0) crossflow line 
at X/L * 0.988 on the floor, (k) la .itudial line at Y/D a 0.683 on the side wall. 



Chapter 6 

CONCLUSIONS AND RECOMMENDATIONS 


Computational simulations of deep cavities in supersonic regimes have been 
performed. Two three-dimensional numerical models have been used to analyze the 
UTisteady flow characteristics of cavity flows. The validity of the codes have been tested 
by analyzing turbulent flows past cavities. Computational simulations of the self 
induced oscillatory flows have been generated through time accurate solutions of the 
Reynolds averaged full Navier-Stokes equations. These governing equations have been 
solved by using an explicit, finite-difference method, and an implicit, finite-volume 
method. The Reynolds stresses have been modeled using the Baldwin Lomax algebraic 
model with certain modifications. Time averaged and instantaneous results have been 
obtained, and quantitative comparisons with experimental data have been made in terms 
of mean static pressure and frequency spectra within the cavity Computational 
simulations of cavity flows are important, because parametric studies on three- 
dimensional cavity flows can be conducted with relative ease and the time dependent 
properties as well as time averaged values can be obtained. In addition, the flow 
structure within the cavity can be visualized computationally. These computational 
capabilities are meant to complement the experimental work in obtaining a more 
complete understanding of the flow features within the cavity. 

Comparing the flow structures of Case 1 and Case 2, we can see the effect of 
yaw. When the yaw angle is zero, the length scale which determines the flow structure 
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!• length-tO i>spth ratio. However, when the yaw angle ie nonzero, there ie no longer 
a tingle length scale for the flow. However, the general structure for L/l) < 3 and 
W/D < 6, Is still an open cavity flow. The dominant directions and inclinations of 
the vortices are dramatically different. The number of vortices are also very different 
from the zero-yaw flow. 

The vortex shedding from the cavity opening is captured computationally. This 
is caused by the motion of the shear layer as it moves in and out of the cavity. Since 
this shedding is predicted by the numerical solutions from both codes, it is believed to 
be a physical phenomenon. It is also believed that this shedding influences the cavity 
acoustics significantly. 

For the explicit code, the explicit addition of the artificial viscosity has the 
tendency to smear out the instantaneous pressure values in the cavity floor region. 
lherefs»e, a time accurate numerical scheme that is robust and capable of yielding good 
results without the addition of artificial viscosity, results in better time «t* pon dont 
results. The implicit upwind scheme, which is naturally dissipative, produces better 
time dependent data. In addition, clustering the grids near the points where the instan- 
taneous pressures are picked, also yields better results. 

Three-dimensional calculations are very expensive, because it takes a long 
computational time to get the results. Two-dimensional calculations of the cavity 
centerplane can be used to develop the cavity flow. When the cavity flow is ftilly 
developed (about two characteristic time), then the solution data for the centerplane can 
be treated as the new initial conditions for the three dimensional cavity. In Case 1, 
such a treatment can save 86% of CPU times and in Case 3 and 4 this treatment can 
save 80% of CPU seconds to develop the cavity flows. 

The acoustic phenomena is a very important topic of cavity flows. The acoustic 
waves generate an additional 20 dB sound pressure level jump in the harmonic frequen- 
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cies. Although wo attempted the acoustic analyses of our results for Cases 2*4, they 
have not been completed. An extension of this work should complete this analysis. 

Quantitative comparison with experimental data have been made in terms of mean 
static pressure and acoustic frequency spectra within the cavity. While most of this 
comparison is favorable, the numerical solution of Case 1 appears to underpredict the 
amplitude of th >«rmonic frequencies. This may be attributed to "too much viscosity" 
due to numerical damping and the eddy viscosity, that alters the pressure oscillations. 
A further study using improved numerical damping and eddy viscosity is needed. 

It is also recommended that the temperature distributions, which are computed 
herein, are evaluated to show the heat transfer characteristics in a cavity. This is 
particularly important for high speed flows. 

It is realized that a better comparison of the two schemes can be obtained if both 
codes were run on the same computers. A further study should be conducted with this 


consideration. 
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